Bulk Viscosity and Relaxation Time of Causal Dissipative 
Relativistic Fluid Dynamics 

Xu-Guang Huang 1,2 , Takeshi Kodama 3 , Tomoi Koide 1 , and Dirk H. Rischke 1 ' 2 
1 Frankfurt Institute for Advanced Studies, 
D- 6 04 3 8 Frankfurt am Main, Germany 
2 Institut fur Theoretische Physik, J. W. Goethe- Universitat, 

D- 6 04 3 8 Frankfurt am Main, Germany 
3 Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
C. P. 68528, 21945-970, Rio de Janeiro, Brasil 

Abstract 

The microscopic formulae of the bulk viscosity £ and the corresponding relaxation time rn in 
causal dissipative relativistic fluid dynamics are derived by using the projection operator method. 
In applying these formulae to the pionic fluid, we find that the renormalizable energy-momentum 
tensor should be employed to obtain consistent results. In the leading order approximation in the 
chiral perturbation theory, the relaxation time is enhanced near the QCD phase transition and rn 
and C are related as rn = C/[^{(l/3-c 2 )(e + J P)-2(e-3P)/9}], where e, P and c s are the energy 
density, pressure and velocity of sound, respectively. The predicted £ and rn should satisfy the 
so-called causality condition. We compare our result with the results of the kinetic calculation by 
Israel and Stewart and the string theory, and confirm that all the three approaches are consistent 
with the causality condition. 
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I. INTRODUCTION 



Relativistic hydrodynamics is an important tool to describe high energy flow phenom- 
ena in different areas, such as heavy-ion collisions, relativistic astrophysics and cosmology, 
although its theoretical foundation has not yet been fully established, in particular, when 
dissipative processes are involved [l| . The simplest formulation is a relativistic covariant ex- 
tension of the non-relativistic Navier-Stokes equation. It is, however, known that relativistic 



Mi. 



Navier-Stokes theory suffers the problem of relativistic acausality and instability 

After the seminal work by Israel and Stewart [s| , so far many different approaches of the 
relativistic hydrodynamics which are consistent with causality have been proposed 0, 6 - 12 1 . 



In the fo 



(CDR) [l]. The crucial difference of any CDR theories from the Navier-Stokes theory can be 
characterized by the introduction of finite relaxation times in the definitions of irreversible 
currents. This aspect is somehow overlooked but has an important consequence. That is, 
any relativistically causal fluids will be non-Newtonian in the sense that irreversible currents 
are no longer simply proportional to the corresponding thermodynamic forces while their 
Newtonian counterparts are. 

In hydrodynamics, all transport coefficients such as the shear viscosity and the bulk 
viscosity are inputs and should be determined from a microscopic theory. In the case of 
the Navier-Stokes theory, the coefficients are usually calculated by using two different ap- 
proaches. One is the kinetic approach based mainly on the Boltzmann equation, and the 
other is the microscopic approach using the Green-Kubo-Nakano (GKN) formula. 

Strictly speaking, the kinetic approach is applicable only to rarefied gas and not reliable 
in practice to calculate the transport coefficients for finite density systems. For exam ple, 
the density expansion of the shear viscosity 77 (p) in three dimensional space is given by 13 1 



lowing, we call these theories as the causal dissipative relativistic hydrodynamics 



V(p) = V0 + V1P + V2P 2 In pH • (1) 

What we can calculate from the Boltzmann equation is only the first term 770- This is because 
the Boltzmann equation is the lowest order approximation of the BBGKY hierarchy equation 
and does not contain the information of multiple collisions which violates the important 
assumption of the Boltzmann equation, that is, the molecular chaos. To calculate the 
further expansion coefficients 771 and 772 in the kinetic approach, we have to use, for example, 
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the Bogoliubov-Choh-Uhlenbeck equation which is a generalized Boltzmann equation 13 ) . 
Unfortunately, systematic generalization of the relativistic Boltzmann equation is not yet 
known. 

On the other hand, the GKN formula does not have such a limitation with respect to 
the finiteness of density as far as we know because it is derived from the microscopic theory, 



quantum field theory 



14l | . However, for a relativistic fluid, we have to be careful because the 



Newtonian property of a fluid is assumed to derive the GKN formula. Thus new formulism 
is needed to calculate the transport coefficients of CDR. 

Transport phenomena such as viscosities, diffusion and heat conduction in hydrodynam- 
ics are related rather to the changes in a sort of boundary conditions than to the responses of 
the system to an external mechanical perturbation. For this reason, the required formulation 
should be different from the ordinary linear response theory. This fact is already emphasized 
by Kubo [is| - As a matter of fact, the well-known expressions of transport coefficients of 
relativistic hydrodynamics are obtained by using the non-equilibrium statistical operator 
method by Zubarev [xsl, • So far, several different approaches have been proposed to cal- 
culate the transport coefficients : indirect Kubo method, Langevin-Fokker- Planck method, 
regression hypothesis based method, local equilibrium approach, external reservoir method, 
and prediction theory, and so on. See Ref. 18| for details. 

Recently, the authors of the present paper proposed a new microscopic formula to calcu- 
late the shear viscosity and the corresponding relaxation time of CDR using the projection 

operator method, which belongs to the Langevin-Fokker- Planck method in the classification 
□ 

mentioned above pj(. Our formula is expressed in terms of the time correlation functions 
of microscopic currents and is a natural extension of the GKN formula. We showed that for 
the shear viscosity, it reproduces the GKN results in the leading order. When applied to a 
Navier-Stokes fluid and non-relativistic diffusion coefficients, our approach reproduces the 



well-known results, as is discussed in Ref. 



20 



2l|. 



The purpose of this paper is to derive the microscopic formulae of the bulk viscosity 
( and the corresponding relaxation time rn in the framework of the projection operator 



method 



22] . We apply the result to the pionic fluid and calculate in the leading order 



approximation in the chiral perturbation theory with the renormalizable energy-momentum 
tensor. We find that the relaxation time rn is enhanced around the temperature near the 
QCD phase transition. 
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We further discuss the differences among our formalism, the kinetic theory [5|, |23H25| 



and the string theory 



B, BE 



26]. In a CDR, the values of ( and rn should satisfy the causality 



condition, which is derived by employing that the propagating speed of a physical signal 
should not exceed the speed of light. If this condition is violated, relativistic fluids becomes 
unstable {3, 4]. We confirm that the values of transport coefficients obtained by all the three 
different approaches are consistent with the causality condition. 

It should be noted that the bulk viscosity and the corresponding relaxation time are 
important not only for heavy-ion collision physics, but also for astrophysics, for example, 
the stability windows in parameter space of rotating compact stars [27]] , driving inflation in 



early universe and the associating entropy production 

This paper is organized as follows. First we give a brief review of the projection operator 
method in Sec. [TT] for the sake of the later convenience. In Sec. IHIt we apply the method 
and derive expressions for the bulk viscosity and the relaxation time. We calculate the 
formulae in the leading order approximation and discuss the consistency with the causality 
condition in Sec. IIV1 Possible other generalization of our result is discussed in Sec. |VJ The 
summary of our results concludes this work in Sec. IVII Through this paper, we will use 
metric g = diag(+, — , — , — ) and the natural units H — c — fee = 1- 



II. PROJECTION OPERATOR METHOD 



For the later convenience, let us briefly review the projection operator method 29N31|. It 



should be emphasized that the projection operator method was firstly proposed by Nakajima 



32J, although it is often referred as the Mori-Zwanzig formalism due to the extensive use 



and developments done by these authors 



33, 



34] 



Many dynamical variables of practical interest, such as the conserved quantities, usually 
vary slowly in time comparing to other microscopic quantities. We call them gross variables. 
1 In order to discuss the dynamics of these slowly varying relevant variables, we need to 
introduce a coarse-graining procedure to smoothing out the microscopic dynamics. The 
projection operator method provides a systematic way to extract the information of the 



1 When a phase transition is present, in principle, we should consider also the corresponding order pa- 
rameters and soft modes as the candidates of the gross (hydrodynamic) variables. We do not, however, 
discuss such a case in this work 
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relevant coarse-grained dynamics from underlining microscopic theories. 

In the case of a quantum system, the full microscopic dynamics is described by the 
Heisenberg equation of motion, 

d t O(t) =i[H,0(t)]=iLO(t) } (2) 

where O is an arbitrary operator and H is the Hamiltonian. For simplicity, here we have 
assumed that H is independent of time. See Refs. 0, Q for the case of a system with 
time-dependent Hamiltonian. The second equality defines the Liouville operator L. In order 
to project out the irrelevant information associated with variables of microscopic (short) 
time scales, we introduce a time-independent projection operator P and its complementary 
operator Q — 1 — P, which satisfy the following general properties, 

P 2 = p, PQ = QP = 0. (3) 



Wit 



34 



i the help of these operators, the Heisenberg equation of motion can be re-expressed as 



36 



37| 



J^O(t) = e iLt PiLO{0) + / dre lL(t ~ T) PiLQe iLQT iLO(0) 
ot J 

+Qe lLQt iLO(0). (4) 

This is called the time-convolution (TC) equation and its r. h. s. is composed of three 
distinct parts. The first term is called the streaming term and usually corresponds to 
collective oscillations such as plasma wave and spin wave. The second term is called the 
memory term which turns into the dissipation term after a coarse-graining procedure. The 
third term is identified with the noise term after implementing coarse-graining of time, as 
we will see later. Thus this equation is considered as a generalized Langevin equation. As 
a matter of fact, the memory term and the noise term are related through the fluctuation- 



dissipation theorem of second kind 



M 



Note that the TC equation is very general and 
still equivalent to the Heisenberg equation of motion. 

The choice of the most appropriate projection operator depends on the specific properties 
of a given system and also on the coarse-graining procedure which we wish to introduce. 
If we choose the projection operator so as to extract all the relevant gross (in our case, 
hydrodynamic) variables, we can, in principle, derive hydrodynamic equations from the TC 
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equation. For this purpose, we use the Mori projection operator 34J. Let 



A 



An 



\ An J 



(5) 



be a n- dimensional vector formed by n time-independent operators corresponding to the 
gross variables {A^ , where the notation A is used to distinguish the Schrodinger operator 
from its Heisenberg form, A = A (t). We choose A = A (0), so that A (t) = e lLt A. 
Then the time-independent Mori projection operator P is defined as 

n 

i=i, 

where O is an arbitrary operator, and the coefficient q is given by 



34] 



(6) 



J2(0,A])-(A,A^ 



3 1 ' 



(7) 



i=i 



Here (X, Y) denotes the inner product of two arbitrary operators X and Y (see below), and 
(A, At)^ denotes ji element of the inverse matrix of (A, A^), i.e., 



5](A ) A%. 1 .(^ J 4) = S itk . 



(8) 



In this way, we expect that the relevant part of an arbitrary operator is expressed as a 
function of the gross variables by operating this projection operator. In order to follow the 
dynamics of the gross variables in time, we would need the time- dependent projections 



but for the present purpose of calculating the transport coefficients, the time independent 
projection is sufficient. 



We are still le 
following Ref. 



t with the freedom to choose the definition of the inner product. Here, 



34J, we use Kubo's canonical correlation, 

{X,Y)= / ^Tr[ Peq e™Xe-> H Y], 
Jo P 



(9) 



where p eq = e ^ H /Tr[e /3H ] with j3 being the inverse of temperature. One can see that, if it 
is a classical system, Kubo's canonical correlation is reduced to the usual classical thermal 
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expectation value. Thus Kubo's canonical correlation is the quantum generalization of the 
classical expectation values. It is easy to confirm that 

(iLX,Y) = -(X,iLY). (10) 

Finally, the TC equation for the gross variable A (t) in the Heisenberg picture can be 
expressed as 

-A(t) = *A A(t) - J drS(r)A(t - r) + £(f), (11) 

where A and H are operators of (n x n)matrices and £ is a n-vector whose elements are 
given by 

•Ay = £(*L^ ) 4)(A ) A%\ (12) 

S ..( t ) = -e(t)^(iLg e * £ «*iLA i> 4)(A > At)^ 1 , (13) 

= Qe tLQt tLA. (14) 

If the set of n-gross variables {A{\ is appropriately chosen so as to extract all the dy- 
namics associated with the slow hydrodynamic time scale, we expect that the dynamical 
variation time-scale of the last term of Eq.f rTTT) should be very small compared to the 
hydrodynamic time scale, because the projection operator Q projects out components only 
orthogonal to {Ai}. For this reason the term £j(t) is called the noise term. 



III. GENERAL FORMULAE FOR BULK VISCOSITY AND RELAXATION TIME 

In this section, we derive the microscopic formulae for the bulk viscosity £ and the cor- 
responding relaxation time rn- Our strategy is as follows. We derive the evolution equation 
of the bulk viscous pressure II from the TC equation, and compare the derived microscopic 
equation with the phenomenological one to extract the microscopic formulae. For our pur- 
pose of obtaining the transport coefficient, it is sufficient to consider small deviation from 
the stationary back-ground-fluid in thermal equilibrium. 

The phenomenological equation of the bulk viscous pressure II in a CDR is given by 0, Q] 

m^n + n = -C9X. ( 15 ) 

where (, rn and u^ 1 are the bulk viscosity, the relaxation time and the fluid velocity, respec- 
tively. The first term on the l.h.s. represents the retardation effect of II which is necessary 
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to satisfy relativistic causality. For rn = 0, Eq. (Tl5|) is reduced to the usual Navier-Stokes 
constructive equation. In general, as is predicted from the kinetic theory, it is possible to 
introduce more non-linear terms in Eq. (1151) but, as mentioned above, here we discuss only 
the lowest order equation consistent with a CDR. 

In order to avoid the possible influence from the shear stress tensor, we consider a pertur- 
bation in infinite fluid in thermal equilibrium having a planar symmetry in (y, z) plane. All 
the quantities associated with the perturbed fluid dynamics vary spatially only along the x 
direction. In this case, the fluid velocity points to the x direction (if one wants to discuss 
shear viscosity only, one can choose the fluid velocit y to point to the x direction but vary 



spatially along the y direction, as is done in Ref. 19|, |29|). Then, the equation of continuity 



of the energy-momentum tensor T^ u in momentum space is given by 

d t T x0 (k x ,t) = -ik x T xx (k x ,t), (16) 

where k x denotes the x-component of the momentum vector k. On the other hand, Eq. (fT5|) 
is simplified as 

T U d t U{k x , t) + U(k x , t) = -(ik x u x {k x , t). (17) 

We use this equation as the definition of ( and r n . Here Tl(k x , t), T x0 (k x , t) and T xx (k x ,t) 
are the Fourier transforms of U(x,t),T x0 (x,t) and T xx (x,t), respectively. 

In order to obtain the microscopic expressions of ( and rn, we derive the equation for 
II from the TC equation. For this purpose, we have to choose appropriate gross variables 
included in Eq. ([17]) to define the projection operator. Among T x0 , T xx , IT and u x , the bulk 
viscous pressure II is the deviation from the equilibrium pressure in the diagonal components 
of the energy-momentum tensor. Thus we use the following operator representation, 

m x ) = -l^(k x ) + 1 -(TH(k x )) eg , (is) 

where ( • • • ) eq represents the equilibrium expectation value. Furthermore, as we will discuss 
later, u x can be regarded as linearly dependent of T 0x in the lowest order of the perturbation. 
Then, Eq. ( fl71) contains basically two independent gross variables, 



A(k x )= / , (19) 




where the bar notation, for example II, refers to the operator value of II (t) at t — 0. 
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According to Eq. ([6]), the projection operator P is now defined by 

(f f °*(-fc*)) 1 j + (n(fc-), n(-^)) 1 J ' 1 j 



Substituting it into Eq. ( TTTj) . we obtain the following two equations, 

d t T 0x {k x ,t) = -ik x U(k x ,t), (21) 

3tII(fcV) = -ik x R k ,T 0x (k x J t) - [ drE 22 (k x ,r)U(k x ,t-r) + C(k x ,t), (22) 

Jo 

where 

fc * (f te (^),f te (-^))' 1 ' 

and we used 

(<LA,At)(A,At)- 1 = f ° ~ lkX 

\ -ik x R k * 

Note that we consider homogeneous energy density and pressure. Then the first equation 
is nothing but just Eq. lfloT) . that is, the equation of continuity, and the second equation 
describes the non-trivial evolution of II. 

The exact expression of the memory function S22 is given in Ref. 



381 ] . However, if we 



are interested in the expression in the low k x limit, we can calculate it more easily. From 
Eq. (J22D, the evolution of (U(k x ,t),U(-k x )) is given by 

<9 t (n(^,t),n(-F)) = -ik x R^{T 0x {k x ,t),Il(-k x ))- [ dTE 22 (k x ,t-r)(U(k x ,r),tl(-k x )). 

Jo 

(25) 

Here we used {£,{k x , t),Yl(—k x )) = 0, which is calculated from the definition of £(k x , t). Then 
the Laplace transform of the memory function at low k x is 

^■^^Sr- < 26 > 

where 



X h (k x ,s) = dte- st \±, , ' , -, K , \ \ (27) 







From the final value theorem of the Laplace transform, we can show that 

v Lnx \ (n(P,t),H(-^)) 
hm sX L (k x ,s) = hm \- - - ' - , -{ L 

(U(k x , 00)) eq (T[(-k x )) eq 



(U(k x ),U(-k x )) 
0. (28) 



Here we used the mixing property of the ergodic theory. Finally, the memory function in 
low k x and s limit is given by 

Z&(K b >s) = ^t£ V- ( 29 ) 

22K ' X L (k x ,s) v ' 



To extract the phenomenological equation (IT7I) . we have to violate the time reversal 
symmetry. For this purpose, we implement the coarse-graining of time. Let us introduce a 
macroscopic time scale tm as follows 

t m = et, (30) 

where e is a scale parameter and less than one. Then the time-convolution integral is 
expressed as 

/ drE 22 (k x ,T)U(k x ,r M /e-r). (31) 
Jo 

When the microscopic and macroscopic time scales are clearly separated, we can take the 
the vanishing e limit. Then the integral is given by 

/•Tivf/e ;>oo 

lim / drE 22 (k x J T)U(k x J T M /e-r)= dr~ 22 (k x , r)U(k x , t) . (32) 
€ ^°Jo Jo 

We call this coarse-grainings the time-convolutionless (TCL) approximation. Note that this 

approximation is very similar to the so-called Markov approximation. In the present case, 

however, there is still the memory effect for II even after the TCL approximation and we 

cannot call it the Markov approximation. With this approximation, Eq. (1221) is expressed 

as 

POO 

d t U(k x ,t) w -ik x R k ,T Qx {k x ,t)- / drZ 22 (k x ,T) U(k x ,t) 

Jo 

= -tk x R k *T 0x (k x ,t) - -^L-^(k x ,t) 

« -tk x R k ,(e + P)u x (k x ,t) - -J—Tl(k x ,t). (33) 
Here, the noise term is neglected. The function r n (k x ) in the second line is defined by 

r n {k x ) =X L (k x ,s = 0). (34) 
From the second line to the third line, we used the following replacement 

T 0x (k x , t)~{e + P)u x (k x , t), (35) 
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which comes from the expression of the phenomenological energy- momentum tensor, 

T 0x (x, t) = [e + P + U(x, t)]u x (x, t), (36) 

and is justified near the local rest frame. Because we defined the projection operator with 
IT and u x by neglecting non-linear terms, we cannot predict the coefficients of non-linear 
terms correctly in the present calculation. Thus, for the sake of consistency, we neglect the 
non-linear term U(x , t)u x (x , t) in Eq. ( 136]) . The validity of the TCL approximation and the 
general comment for the derivation of the non-linear term are discussed in Sec. [V] 
By comparison Eq. (|33|) with Eq. ( ITT)) , we obtain the following correspondences, 



r n = lim X L {k,s), (37) 
( = ( £ + P)R lim X L (k,s). (38) 

For the sake of convenience, we express these expressions in terms of the retarded Green 
functions. Note that the correlation function X L (k, s) can be re-expressed as 

X L (k,s) 

-i -l 



1 

1 Jo 



OO /• fOO 

j3„ st— ikx 



dt / cfxe- s *- Jk - x / dr(n(x,r)n(0)) ret 



d 3 xie- <k - xi (n(xi),n(o)) 



(39) 



where the retarded Green function is defined by 



(n(x,t)n(xi,r)) rrf = -ze(t-r)Tr{p ef/ [n(x,t),n( Xl ,r)]} 

= /_^^Gg(c,k)e-('-)e-^™). (40) 

See Appendix |A] for details. 

Then, finally, the formulae of ( and rn are given by 

C Cgkn 



f3{e + P) /3 2 /c/ 3 x(f te (x),f to (0)) 
■ lim lim dG^(u, k)/du 



13 lim lim G% x (lo, k) 



GKN 



(3 /3 2 /d 3 x(n(x),n(o)) 



■ lim lim <9Gff(u;, k)/du 
"P lim lim gWcu, k) 

k->0 w^O 
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(41) 
(42) 
(43) 
(44) 



Here we have introduced the usual expression of the bulk viscosity for the relativistic Navier- 
Stokes fluid in the GKN formula (obtained with the Zubarev method), 



GKN 



/poo poo 
d 3 x J dt J rfr(n(x,r)n(0)) re4 



u->0k->0 duo 



and one more retarded Green function, 



/oo 
rfx(T 0x (x, t)r te (0)) rei e- iw *e <k,c . (46) 
-oo 

Equations (|4ip -(l44 p are our main results. The bulk viscosity and it relaxation time are 
expressed by the ratios of Green's functions and different orderings of limits. 



IV. APPLICATIONS TO HOT PIONIC FLUID 

As an application of our microscopic formulae (142]) and ( |44|) . we will calculate the bulk 
viscosity Q and relaxation time rn for hot pion fluid in confined phase within an effective 
model. 

Let <p be the scalar field for pions (we simply use a real scalar field to present pions since 
the charge does not affect the results). The usual definition of the energy- momentum tensor 
of this field is 

= d»(j)d"<f) - g^ u C, (47) 

where £ is a Lagrangian density. In this case, the bulk viscous pressure ffl8|) for non- 
interacting case, would become 

LT(x,t) = i[(90(x,t)) 2 - 2MVW)] - i((<9<Kx,t)) 2 - 2M 2 2 (x,t)> e „ (48) 

where M is the mass of pion. However, these expressions are not adequate for our purpose. 
First, note that the above bulk viscous pressure does not vanish even in the massless limit 
M = 0, which does not reflect the conformal property of the Lagrangian in this limit. Fur- 
thermore, the energy- momentum tensor (T47j) is not renormalizable, i.e., its matrix elements 
depends directly on the cut-off of the renormalized perturbation theory as is discussed in 
Ref. |40| . Thus we introduce the renormalizable energy- momentum tensor 6^ following 
Ref. (40j as, 

^(x, t) = T^(x, t) - - ^9 2 )0 2 (x, t). (49) 
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Then the corresponding bulk viscous pressure for the non-interacting case is given by 



1 



M 2 



n(x,t) = --0£(x,t) + -m*,t)U = -— (0 2 (x,t) 



(M)> 



e?7 ) 



(50) 



which recovers the conformal nature of the system in the vanishing limit of M. Note here 
that, for fermion field and gauge field, the usual definition of energy- momentum tensor is 
already renormalizable and do not need any re-definition of the energy-momentum tensor. 

However, because of the reason which will be discussed in the end of this section, this is 
still not the definition of the bulk viscous pressure which is used in the following calculation. 
We recall that the behavior of the retarded Green function in the low momentum limit 
is not changed by adding an additional term which is proportional to the energy density in 
the definition of the bulk viscous pressure. Finally, we added an additional term which is 
proportional to the energy density to define the bulk viscous pressure instead of Eq. (fl8l) . 



n(x,t) 



i 



^(x,t)4 
M 2 



3 ^'' J ' \3 



# 00 (x,t) 



i-c 2 

3 s 



3 ^ ' \3 
^(^V(x,t) 



— c: 



# 00 (x,t) 



eq 



i-c 2 
3 s 



9™(x,t)) (51) 



ey 



where c s is the velocit y o f sound. This is the same definition of the bulk viscous pressure 



discussed in Refs. 



16 



4l|. One can easily see that this bulk viscous pressure still vanishes 



in the massless limit. 

As the lightest particles, pions dominate the transport properties of QCD in hadronic 
phase. From Eqs. ( jUJ) and (l4"3|) . once the leading order result of (gkn is obtained, the cor- 
responding leading order results for ( and Tn are obtained by substituting the denominators 
on r.h.s by their non-interacting counterparts. A straightforward calculation leads to 



lim lim G% x {us, k) = — e — P, 



HmlimGg(u;,k) 



1 A / „x e-3P 



(52) 
(53) 



where E p = ^/p 2 + M 2 . The energy density and pressure of the free pion gas are, respec- 
tively, given by 



N, 



V 



3^ ^ R 



(54) 
(55) 
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where f(x) is the Bose-Einstein distribution function \/{eP x — 1) and the prefactor N n = 3 
counts the degeneracy of ir + , 7r~ and 7r°. 

The r. h. s of Eq. (l53p has a term which contains ultraviolet divergent vacuum term. In 
order to obtain the finite result, we renormalized this vacuum term. However, it should be 
noted that the renormalization of this divergence is not trivial as is discussed in Ref. 4^ |. 

In short, in the leading order approximation, the bulk viscosity and the relaxation time 



are given by 



22) 



GKN 



(3{e + P) 0(e + P) 



Cgkn 



P) 



3P 
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The first equation shows that the bulk viscosity ( is reduced to the GKN bul 
Cgkn-i similarly to the case of the shear viscosity in the leading order calculation 



(56) 
(57) 

£ viscosity 



19]. There 



ravior here. 



44| , we plot 



already exist several calculations for Cgkn [43|. Thus we will not discuss its be 
By adopting the result of (gkn calculated in the chiral perturbation theory 
the temperature dependence of the dimensionless ratio Tn//3 in the hadron phase in Fig. 
[TJ For the sake of comparison, the relaxation time calculated from the Boltzmann equation 
with Grad's moment method is shown by the dashed line 45[. The order of the magnitude 
is same as the relaxation time of the shear viscosity, r n //3, which is shown in Fig. 1 in 



Re 



. 19[. However, the temperature dependence of 7~n//3 is non-trivial. As is shown in Ref. 
191 ]. r n /P is a monotonically decreasing function of temperature in the hadronic phase. On 
the other case, 7~n//3, which is a decreasing function at low temperature, starts to increase 
around T = 100 MeV and shows maximum near the QCD phase transition. This comes 
from the enhancement of (gkn caused by the trace anomaly. See Ref. for details. 



Now we compare our result with the resu 
approximation [5 



23] and the string theory 



ts from Grad's method with the 14 moment 
. For this purpose, it is convenient to consider 



the (-tjj ratio, because this quantity is independent of the choice of the collision term in the 
Boltzmann equation. In our leading order result of pion, this ratio is given by 



c 



R 



(l-c 2 s ) (e + P)-l(e-3P) 
(e + PA 



In the result of the string theory, this ratio is given by 



26) 



(1/3 



Tn(e + P) 2 -In 2 



(58) 



(59) 
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FIG. 1: The temperature dependence of the relaxation time of the bulk viscous pressure rn//3 of 
the hadron phase. The solid line represents the result of the leading order approximation in the 
chiral perturbation theory. For the sake of comparison, the result of the 14 moment approximation 



from Ref. 



451 ] is shown by the dashed line. 



In the 14 moment approximation, this ratio is calculated by using a function /3q which is 
denned by Eq. (7.8c) of Ref. I|. 

The temperature dependence of the (-Tu ratio for pion is shown in Fig. [2J The solid, 
dashed and dot-dashed lines represent our microscopic formula fl58l) , the 14 moment approx- 
imation, and the string theory, respectively. The string theory predicts the largest value of 
the ratio, meanwhile the 14 moment approximation estimates the smallest value. At high 
temperature, the three lines are monotonically decreasing functions of temperature. The 
qualitative difference is observed at low temperate. The ratio of our formula and the string 
theory are still finite at T = but that of Grad's method vanishes. The meaning of this dif- 
ference is related to quantum fluctuation which is not included in the Boltzmann equation. 
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FIG. 2: The temperature dependence of the C" T n ratio. The dot-dashed, solid and dashed lines 
represents the results of the 14 moment approximation, our formula (|58p and the string theory, 
respectively. 

This result is reported in another paper comparing our result with a new kinetic calculation 



based on the Boltzmann equation 22 



It is worth for mentioning that this ratio is closely related to the propagation velocity of 
signals in CDR, which is given by {3, 7] 



c 



(60) 



iere, we neglected the contribution from the shear stress tensor. As is discussed in Refs. 
3|, |4|, for hydrodynamics being causal and stable, this group velocity should be smaller than 
the speed of light. Thus this ratio should satisfy the following constraint, 

c 



Me + P) 



< 1-c 



(61) 
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This is the so-called causality condition {3), 0, One can easily see from Fig. [2] that all the 
three calculations satisfy the causality condition. 

As was mentioned before, we have used Eq. ( I51|) instead of Eq.f l50p as the definition of 
the bulk viscous pressure. Then we pointed out that both definitions should give the same 
result. However, this is not trivial when there is a UV divergence. As a matter of fact, 
the C/ r n ratio calculate with Eq. ( )50|) is different from the solid line in Fig. [2] when we 
adapt the simple subtraction of the vacuum term as renormalization. In this calculation, we 
believe the calculation with Eq. (l5Tj) is more reliable than that with Eq.( l50|) because of the 
following three reasons. 1) When we use Eq. ( I5U|) and employ the simple renormalization, the 
calculated C/ r n violate the causality condition. 2) As is pointed out in Ref. 4l|, if Eq. flBTJj) 
is used, the perturbative calculations collapses because of divergence. This problem is solved 
by using Eq. (I51[) . 3) We can show that the C/ r n calculated with Eq . floT]) is consistent with 



the result from a new kinetic calculation, as is discussed in Ref. 



22]. 



V. OTHER POSSIBLE GENERALIZATION 



A. Another approximation to the memory function 



In the derivation of Eq. ( l33l) . we replaced the time-convolution integral of the memory 
term with the time-convolutionless integral by assuming that the macroscopic time scale is 
clearly separated from the microscopic one. On the other hand, this coarse-grainings may 
be formulated as a following expansion of the memory term, 



drE 22 (k x ,T)Il(k x ,t - t) = / drE 22 (k x ,T) 
Jo 



Tl(k XJ t)-T^Tl(k x ,t) + 



(62) 



When the higher order correction becomes important, the evolution equation of the bulk 
viscous pressure is modified as 



drZ 22 (k x ,T) 
— > f n d t U(k x ,t) + U(k x ,t) = -ik x (u x (k x ,t). 



U{k x ,t)-T^U(k x ,t) 



(63) 



Here, we replaced the upper limit of the integral by 00, assuming that the dominant contri- 
bution of the memory function still comes from r = 0. The expression of ( is not changed 

by this correction. However the relaxation time is modified by fn which is given by 

8E L {k, s) 



Tn 1 + lim 

s,k->0 



ds 



(64) 
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This can be expressed in terms of the retarded Green function as follows, 

( lim d 2 G^(co,k)/duj 2 )( lim G*{u,k)) 
Tn = Tn 2( hm 9Gg(a;,k)/^) 2 " (65) 

For example, let us assume an exponential form for the memory function, 

S 22 (0,t) = ATe- rt , (66) 

where A and T are parameters. Then we obtain that 

Tn = r n (l - . (67) 

The result with the TCL approximation is reproduced in the infinite T limit. 

However, because of the following reasons, we consider that the rn is more reliable than 
fn- Let us consider the following equation, 

d t J{t) = F{t)+ [ drZ{t-T)J{T). (68) 
Jo 

From the initial and final value theorem of the Laplace transform, we can calculate the 
initial and final values of J(t) as follows, 

lim J(t) = lim sJ L (s) = lim g = °) ~ ^ L ( g ) = j< t = Q) _ F ^( 00 ) ) ( 69 ) 

t->0 s— >oo s—>-oo s — C (s) 

lim J(t) = limsJ L (s) = 0. (70) 

Here we assumed that S L (0) is finite because of the existence of the finite relaxation time. 
Next, we approximate the time-convolution integral of the above equation by using the 
Taylor expansion, 

d t J(t) = F(t) + AJ(t) + Bd t J{t). (71) 
Then the initial and final values are 

= J( f = 0)-BA f = 0)-^(co) 
t^o w 1 — B 

lim J(t) = limsJ L (s) = 0. (73) 

t— >oo s— >0 



In order to reproduce Eq. (1691 . we have to set B = 0. 

As another reason, we consider the calculation of the transport coefficients of the diffusion 
equation. As was shown by Kadanoff and Martin j^J and later by one of the present authors 

18 



[47(, the ratio of the diffusion coefficient D and the corresponding relaxation time td is 
exactly determined from a sum rule. See Eqs. (B.10), (B.16) and (B.28) of Ref. [4^] . Then 
we obtain 

D_ = Jd 3 x(J(x),J(0)) . . 

t d /d3x(n(x),n(0))' 1 J 

where n is the conserved number density and dtn + VJ = 0. In the projection operator 
method, this result is reproduced only when the TCL approximation is applied. 

That is, if we regard that the TCL approximation as the lowest order of the Taylor expan- 
sion and the next order correction is considered, we obtain the result which is inconsistent 
with the initial value theorem and the sum rule. 

When we observe with the time scale where the structure of the memory function in Eq. 
(168j) cannot be neglected, it means that there still exists macroscopic degrees of freedom 
in the memory function. In order to finish the program of coarse-graining in the projec- 
tion operator method, we have to re-define the projection operator so as to extract this 
macroscopic degree of freedom. Then the form of hydrodynamics itself is changed and the 
obtained equation is not given by Eq. (|15p anymore. In other words, once we assume 
that the hydrodynamic equation is given by Eq. (I15jl . we must observe with the time scale 
where the structure of the memory function is negligible, and hence, we must use the TCL 
approximation instead of the Taylor expansion. 



As an example of the time dependence of the memory function, see Fig. 4 in Ref. 35]. 
One can see that the memory function has finite values only around t = 0. 

B. Non-linear terms 

In this work, we used the phenomenological equation l fl5l) as the definition of the bulk 
viscosity £ and the relaxation time tu- On the other hand, it is possible to introduce non- 
linear terms to the phenomenological equation. For example, it might be possible to derive 
the hydrodynamic equation for the rarefied gas from the Boltzmann equation, by using 
Grad's moment method with the 14 moment approximation. In particular, the equation of 
the bulk viscous pressure is given by, 



23 



2- 



r n —U + II = -C<9X + ain<9X + + + a^A^a + a^A^daUp, 



(75) 
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where 



= g" u - u**u u , (76) 

^fli/aP _ \ ^^fia^iyp _|_ _ ^^u^a^ ^jtj^ 

In order to obtain the microscopic expressions of these a^s in the projection operator method, 
the definition of the projection operator must be generalized so as to collect all gross variables 
appears in this non-linear equation. This will be put in future work. 

VI. CONCLUDING REMARKS 

In this paper, we derived the microscopic formulae of the bulk viscosity £ and the cor- 
responding relaxation time tjj in causal dissipative relativistic fluid dynamics by using the 
projection operator method. Applying these formulae to the pionic fluid and calculating in 
the leading order approximation in the chiral perturbation theory, we found that the relax- 
ation time is enhanced around the temperature near the QCD phase transition, and there is 
a simple relation rn = (/[/3{(l/3— cf)(e+P) — 2(e — 3P)/9}] between ( and rn- We compared 
our result with the results of Grad's moment method with the 14 moment approximation 
and the string theory by calculating the (-Tn ratio which is independent of the choice of the 
collision term of the Boltzmann equation. The ratio must be smaller than (1 — c 2 s )(e + P) 
to satisfy the causality condition. We confirmed that all the three approaches are consistent 
with the causality condition. Finally, we discussed that the time-convolutionless approxi- 
mation, which is used to derive the transport coefficients, is consistent with exact results 
and we should not consider corrections to this approximation. 

It should be emphasized that Grad's moment method with the 14 moment approximation 
is not the unique method to calculate the transport coefficients of CDR from the Boltzmann 
equation. Recently a new calculation method based on the Boltzmann equation was devel- 
oped and it was found that the calculated transport coefficients are different from those of 
the 14 moment approximation {25)]. These new results are completely consistent with the 



leading order results of our formulae. This result was reported in another paper [22]. 

We thank T. Brauner, G. S. Denied, D. Fernandez-Fraile, J. Noronha, H. Warringa, 
and Z. Xu for useful discussions and comments. This work was (financially) supported by 
FAPERJ, CNPq and the Helmholtz International Center for FAIR within the framework 
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of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich- Okonomischer 
Exzellenz) launched by the State of Hesse. 

Appendix A: Relation between Kubo's canonical correlation and retarded Green's 
function 

In this appendix, we give a relation connecting the Kubo's canonical correlation and the 
retarded Green's function. This relation is used to obtain Eqs. ( l4"2"j) and (|44p . From the 
definition of Kubo's canonical correlation, we can transform as 



(n(t,x),n(t',x')) = 



(3 



Tr [p eq e XH U(t,x)e- XH Il(t' , x' 



)] 





f" dX f°° , /TT/ , , J „, 

r°° l"P rl\ rt 



S ,x')} 



eg 



S ,x')> 



eg 



(Al) 



Furthermore, when t = t', we have 



f°° r/s 

(n(o,x),n(o,x')) = • / -e( s )([n( s ,x'),n(o,x)]} 



eg 



«/ — oo 




(A2) 



Then we obtain 



d 3 x(IT(0,x), 11(0,0)) 




'fi(^k). 



(A3) 
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